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ABSTRACT 

Context. A bright X-ray transient was seen during an XMM-Newton observation in the direction of the Small Magellanic Cloud 
(SMC) in October 2006. 

Aims. The EPIC data allow us to accurately locate the source and to investigate its temporal and spectral behaviour. 
Methods. X-ray spectra covering 0.2-10 keV and pulse profiles in different energy bands were extracted from the EPIC data. 
Results. The detection of 6.85 s pulsations in the EPIC-PN data unambiguously identifies the transient with XTE J0103-728, discov- 
ered as 6.85 s pulsar by RXTE. The X-ray light curve during the XMM-Newton observation shows flaring activity of the source with 
intensity changes by a factor of two within 10 minutes. Modelling of pulse-phase averaged spectra with a simple absorbed power-law 
indicates systematic residuals which can be accounted for by a second emission component. For models implying blackbody emission, 
thermal plasma emission or emission from the accretion disk (disk-blackbody), the latter yields physically sensible parameters. The 
photon index of the power-law of ~0.4 indicates a relatively hard spectrum. The 0.2-10 keV luminosity was 2xl0 37 erg s with a 
contribution of ~3% from the disk-blackbody component. A likely origin for the excess emission is reprocessing of hard X-rays from 
the neutron star by optically thick material near the inner edge of an accretion disk. From a timing analysis we determine the pulse 
period to 6.85401(1) s indicating an average spin-down of ~0.0017 s per year since the discovery of XTE J0103-728 in May 2003. 
Conclusions. The X-ray properties and the identification with a Be star confirm XTE J0103-728 as Be/X-ray binary transient in the 
SMC. 
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1. Introduction 

The transient pulsar XTEJ0103-728 with a period of 
6.8482+0.0007 s was discovered during RXTE observations of 
the Small Magellanic Cloud (SMC). The pulsations were de- 
tected on 2003, Apr. 29, May 7, May 15 a nd May 19 but not 
■ on Apr. 24 and May 28 dCorbet et alj|2003l) . suggesting an X- 
ray outburst which lasted about three to four weeks. This be- 
haviour is characteristic for Be/X-ray binaries undergoing a type 
II outburst which is caused by the ejection of matter from the 
Be star and enhance d accretion onto a compact object, in most 
cases a neutron star dNegueru ela 1998). Shorter outbursts (type 
I) with a duration of a few days are often separated by the or- 
bital period of the binary system. They generally occur close to 
the time of periastron passage of the neutron star when the neu- 
tron star approaches the circu mstellar disc of the Be star (see e.g. 
lOkazaki & Negueruelall200ll) . 

Be/X-ray binaries form the major class of High Mass X-ray 
Binaries (HMXBs). In the remaining systems - the supergiant 
HMXBs - the compact object accretes matter from the fast stel- 
lar wind of an early type O or B supergiant star. In the SMC 
more than 60 Be-HMXBs are known, while only one HMXB is 
established as supergiant system (SMCX-1). Recent reviews of 
th e optical and X-r ay p roperties of these system s can be found 
in lCoe etal.1 d2005l) and lHaberl & Pietschl d2004l) . respectively. 

* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by ESA 
Member states and the USA (NASA) 



The X-ray transient XTE J0103-728 is one of several Be/X- 
ray binary candidates in the SMC which were discovered 
through their pulsations with RXTE. This non-imaging instru- 
ment is able to mo nitor the activity of pulsars in the SMC 
dLavcock et al.ll2005l) but provides only a coarse determination 
of their sky coordinates. The latter hampers the optical identifi- 
cation and final confirmation of the HMXB nature of the pulsars. 
The pulsar XTE J0103-728 was detected in X-ray outburst dur- 
ing an XMM-Newton observation of RX J0103.8-7254, a can- 
didate super-soft X-ray source in the direction of the SMC dis- 
covered by ROS AT. The precise localization in the EPIC images 
allowed us to identify the optical counterpart of XTE J0103-728 
which sh ows optical brightn ess and colours consistent with a Be 
star (see lHaberlet all 120071) . Analysis of MACHO and OGLE 
data of the cou nterpart reveals long-term va riations in the opti- 
cal light curve (Schmidtke & Cowlevl 120071) on time scales of 
620 to 660 days which are not periodic, suggesting that they 
are related to the Be p henomenon and not to the binary orbit 
dMcGowan et al.ll2007l) . 

Here we present a temporal and spectral analysis of the EPIC 
data of XTE J0103-728 obtained in October 2006. 

2. Observations and data reduction 

An X-ray transient was seen a s brightest source in the EPIC field 
of view of an XMM-Newton ( Jan sen et alJ uOOl) observation in 
the dir ection of the SMC, as summarized in Table [T] The EPIC- 
MOS dTurner et alj l200l and EPIC-PN dStriider et alj l200lh 
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Table 1. XMM-Newton EPIC observations of XTE J0103-728. 



Observation 


Pointing direction 


Sat. 


ID 


R.A. 


Dec. 


Rev. 




(12000.0) 




0402000101 


01 03 52.2 


-72 54 28 


1248 


EPIC 1 


Start time 


End time 


Net 


instrument 


2006-10- 


03 (UT) 


exposure 


configuration 






ks 


PN FF thin 


00:31:21 


06:13:06 


17.72 


Ml FFthin 


00:08:39 


06:12:46 


21.32 


M2 FF thin 


00:08:39 


06:12:51 


21.30 



1 FF thin: full frame CCD readout mode with 73 ms frame time for PN 
and 2.6 s for MOS; thin optical blocking filter. 



cameras were operated in imaging mode covering the source 
at an off-axis angle of about 12'. Background flaring activity 
was negligible and we used the full exposure time for our analy- 
sis. For the data processing we used the XMM-Newton Science 
Analysis System (SAS) version 7.1.0 supported by tools from 
the FTOOL package together with XSPEC version 11.3.2p for 
spectral modelling. 

3. Results 

We performed a source detection analysis of the EPIC im- 
ages using standard maximum likelihood procedures from the 
SAS package. After applying a bore-sight correction using a 
background AGN and a foreground star in the field of view 
the position of the brightest source was determined to R.A. 
= 01 h 02 m 53 s .39 and Dec. = -72°44'3476 (J2000.0) with a re- 
maining systematic uncertainty of 1.1" (l<x ). No source at this 
position was seen in ROSAT observations dHaberl et al.l feOOO; 
ISasaki et al.|[2000t) . 

During the XMM-Newton observation the source showed 
strong flaring activity with intensity changes by a factor of ~ 2 
within 10 minutes. This is illustrated in Fig. Q] where the broad 
band EPIC-PN light curve is shown. A Fourier analysis of the 
EPIC-PN data clearly revealed the presence of 6.85 s pulsations 
in the X-ray flux which identifies the bright EPIC source with 
XTE J0103-728. No significant signal of higher harmonics is 
seen. Using a folding technique the pulse period is determined 
to (6.85401 ± lxl0~ 5 ) s (lcr error). This corresponds to an av- 
erage spin-up of ~0.0017 s per year between the RXTE and the 
XMM-Newton observations, 3.4 years apart. 

The precise position al lowed us to identify the optical coun- 
terpart dHaberl et alj2007l) which is included in the UBVR CCD 
Survey of the Magellanic Clouds (lMassey|2002l). the Mag ellanic 
Clouds Photometric Survey (MCPS IZaritskv et al.ll2002l) and in 
the OGLE BVI photometry catalogue ludalski et alJll998l) as 
summarized in Table [2] Optical brightnes s, colors and temporal 
properties are all consistent w ith a Be star ( Schm idfke & Cowlevl 
2007; iMcGowan eTail 120071) . confirming the Be/X-ray binary 
nature of XTE JO 103-728. 

The X-ray light curves were folded in standard EPIC energy 
bands as shown in Fig. [2] Hardness ratios were derived with 
HR1 = (R2-R1)/(R2+R1), HR2 = (R3-R2)/(R3+R2) and HR3 = 
(R4-R3)/(R4+R3) with RN denoting the background-subtracted 
count rate in band N (starting with band 1 at the lowest energies). 
HR1, which is computed from count rates below 1 keV is more 
sensitive to changes in absorption column density while HR2 
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Fig. 1. Broad-band EPIC-PN light curve of XTE JO 103-728. The 
data is background subtracted and binned to 100 s with time 
corresponding to MJD 5401 1.02678. 



and HR3 can be used as indicator for variations in the shape of 
the intrinsic source spectrum up to 4.5 keV. The three hardness 
ratios plotted versus pulse phase are shown in Fig. [3] 

The pulse profile of XTEJ0103-728 is highly structured and 
strongly energy dependent. The pulsed fraction decreases with 
increasing energy up to ~4.5 keV, but is higher again above ~4.5 
keV. Several features in the folded light curves strongly change 
their appearance in the various energy bands. E.g. the intensity 
maximum in the 0.2-0.5 keV band reverses to a minimum at 
high energies above 4.5 keV and the sharp minima seen between 
1.0-2.0 keV are largely smeared out at energies above 4.5 keV. 
HR1 indicates a relatively smooth change in the low-energy part 
of the spectrum, either due to variations in column density or 
in a low-energy spectral component. Variations are also seen in 
HR2 and similarly in HR3 which suggest changes in the overall 
spectral shape. 

We extracted pulse-phase averaged EPIC spectra for PN (sin- 
gle + double pixel events, PATTERN 0-4) and MOS (PATTERN 
0-12) excluding bad CCD pixels and columns (FLAG 0). Due 
to the large off-axis angle we used an elliptical source extraction 
region which was cut in the case of EPIC-PN by the nearby CCD 
border. The fluxes derived from the PN spectrum were system- 
atically lower than those from the MOS, probably due to under- 
correction of the flux losses. Therefore, we report fluxes and lu- 
minosities (see below) derived from the MOS spectra. The three 
EPIC spectra were simultaneously fit with the same model al- 
lowing for a constant normalization factor between the spectra. 
The spectra show an extremely flat distribution in the EPIC en- 
ergy band which indicates a very hard spectrum. We included 
two absorption components in our spectral modelling, account- 
ing for the Galactic foreground absorption (with a fixed hydro- 
gen column density of 6xl0 20 cm' 2 and with elemental abun- 
dances from lWilms et al.l d2000l) ) and the SMC absorption (with 
column density as free parameter in the fit an d with metal abun- 
dance s reduced to 0.2 as typical for the SMC (Russell & Dopita 
|1992|) ). A bremsstrahlung model is not able to reproduce the flat 
spectral distribution while a power-law (PL) yields a formally 
acceptable fit (see Table [3] and Fig. 0). However, the residu- 
als indicate some systematic deviation from a power-law model. 
Using a broken power-law yields a good fit with ^ 2 =1.02 with 
a steeper power-law index below 2 keV than above, which in- 
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Table 2. Optical identification of XTE JO 103-728. 
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Catalogue 


R.A. and Dec. (J2000.0) 


Vmag 


B-V 


U-B V-R 


V-I 


UBVR 


01 h 02 m 53 s :30-72°44'34'.'9 


14.59 


-0.08 


-0.96 -0.90 




MCPS 


01 h 02 m 53 s .39 -72°44'34'.'7 


14.99 


-0.19 


-1.12 


-0.10 


OGLE 


01 h 02 m 53 , :29 -72°44'34'.'8 


14.68 


-0.12 




0.04 




0.0 0.5 1.0 1.5 2.0 

Phase (6.854010 s) 



Fig. 2. Folded EPIC-PN light curves in the standard EPIC en- 
ergy bands. The panels show the pulse profiles for the different 
energies specified in keV. The intensity profiles are background 
subtracted and normalized to the average count rate (in cts s _1 : 
0.103, 0.268, 0.482, 0.530, 0.466, 1.86 from top to bottom). 



dicates excess emission below ~2 keV. For a better description 
of this excess we add different model components to the simple 
power-law. 

Two different kinds of low-energy excess seen in the X- 
ray spectra of HMXBs were reported in the literature. A soft 
emission component from Be-HMXBs is interpreted as re- 
processed X-ray emission and may originat e from different 
places in the binary system dHickox et alj2004l) . The Be-HMXB 
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Fig. 3. Hardness ratios HR1, HR2 and HR3 derived from the 
folded EPIC-PN light curves shown in Fig. [2] 



RXJ0103. 6-7201 (with 1323 s the pulsar with the longest pe- 
riod known in the SMC) showed during one XMM-Newton ob- 
servation a highly absorbed power-law component and a com- 
pletely disentangled soft component. The soft component wa s 
modelled by thermal plasma emission (lHaberl & Pietsch 2005). 
Using a thermal plasma emission component in combination 
with the power-law (PL+MEKAL) yields an improved fit for 
XTEJ0103-728 (Table©, but a temperature of > 1.3 keV much 
higher than typical values found for other sources (0.15 keV for 
RXJ0103 .6-7201). This is inconsistent with a very soft emis- 
sion component present below ~ 1 .3 keV, similar to that observed 
from RXJ0103. 6-7201. A relatively hot component with char- 
acteristic blackbody temperature of il keV and a small emission 
area with radius <0.5 k m was found in two persistent Be pu lsars 
(e.g. RX J0146.9+6121: lLa Palombara & Mereghettill2006l) and 
was interpreted as emission from the hot polar caps of the ac- 
creting neutron star. Adding a blackbody component to the mod- 
elling of XTEJ0103-728 (PL+BB) we obtain a lower tempera- 
ture around 0.25 keV (Table© and a larger emission area (radius 
30.8 km) which is inconsistent with the surface of the neutron 
star. The inferred radius might indicate emission from the in- 
ner accretion disk and therefore, we added a multi-temperature 
blackbody model (diskbb in XSPEC) to the power-law model 
(PL+DB in Table©. Again, the fit is acceptable with a temper- 
ature of the inner disk of 0.36 keV (Table© and an inner disk 
radius of ~ 18.7 km (assuming an average disk inclination angle 
cos 6* of 0.5). The disk-blackbody component contributes 3.0% 
to the total luminosity in the 0.2-10.0 keV band. The EPIC spec- 
tra together with the best fit PL+DB model are shown in Fig|4] In 
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Fig. 4. EPIC spectra of XTE J0103-728. EPIC-PN is shown in 
black and EPIC-MOS in red (Ml) and green (M2) (both grey in 
black and white representation). The histograms show the best- 
fit model: (top) an absorbed power-law and (bottom) an absorbed 
two-component model with multi-temperature disk and power- 
law emission (PL+DB in TableO. 

Table [3] we summarize the characteristic model parameters and 
give observed fluxes and luminosities. 

4. Discussion 

After its discovery with RXTE in April 2003, the 6.85 s pul- 
sar XTEJ0103-728 was detected in X-rays for the first time 
with imaging instruments which allowed us to better localize the 
source and identify the optical counterpart. During the XMM- 
Newton observation in October 2006 the source was seen in out- 
burst with a luminosity of ~2xl0 37 erg s _1 in the 0.2-10.0 keV 
band. During the observation the source strongly varied in inten- 
sity with changes by a factor of two within ten minutes. Pulse 
profiles obtained by folding the data from the EPIC-PN detector 
in different energy bands with 30 phase bins (corresponding to 
228.5 ms per phase bin) are highly structured. The fast intensity 
variations indicate that the bulk of the emission originates close 
to the neutron star surface and suggests a complicated emission 
geometry with narrow beams contributing to the X-ray emission 
at various neutron star spin phases and with energy dependent 
amplitudes. 

Deviations from a power-law model in the phase averaged 
EPIC spectra suggest an additional emission component mainly 
contributing at energies below ~1.5 keV. This component is 
harder than the soft component seen in other Be-HMXBs in the 



SMC, which can be modeled by thermal plas ma emission with 
kT~0 .15 keV (see e.g. RX JO 103 .6-7201; lHaberl & Pieischl 
2005). On the other hand the excess emission in the spectra of 
XTE J0103-728 is softer than the hot blackbody component seen 
in XMM-New ton spectra of two persistent Be -HMXBs in the 
Milky Way dLa Palombara & Mereghetti IIM 120071) which is 
interpreted as emission from hot polar caps of the neutron star. 
When modelled as blackbody emission, the large inferred emis- 
sion area (with radius ~30 km) for XTE J0103-728 is incompat- 
ible with an origin on the neutron star surface. 

In their work on the origin of the soft excess in X-ray pulsars, 
iHickox et al.l (2004) conclude that for luminous X-ray pulsars 
(ilO 38 erg s" 1 ) the soft excess can only be explained by repro- 
cessing of hard X-rays from the neutron star by optically thick 
material. At intermediate luminosity (~10 37 erg s ) also other 
processes such as emission from photo-ionized or collisionally 
heated gas can contribute. The spectra of XTE J0103-728 do not 
show significant line emission suggesting that in this source the 
excess emission is mainly caused by reprocessing in optically 
thick material, most likely located near the inner edge of an ac- 
cretion disk. The relatively low contribution (~3%) to the total 
luminosity and the estimate for the inner disk radius of ~19 km 
inferred from the disk-blackbody model are consistent with such 
a model. 

The luminosity of the soft component seen from 
RXJ0103. 6-7201 is strongly correlated with the total source 
luminosity over at least a factor of ten variation in source in- 
tensity. At a maximum of ~6.4xl0 36 erg s" 1 RX J0103 .6-7201 
is still a factor of ~3 below the luminosity of XTEJ0103-728 
during its Oct ober 2006 ou t burst. This is consistent with the 
conclusions of Hic kox et alj d2004l) . that for lower luminosities 
emission from diffuse, optically thin gas dominates. Which of 
the emission components becomes visible in the X-ray spectra 
of Be-HMXBs certainly depends not only on source luminosity 
but also on geometric effects. 
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Table 3. Spectral fit results. 



Model ' 


SMC N H 


r 


kT 


Flux <2) 


L <3) 


^, 2 /dof 




[10 20 crrT 2 ] 




[keV] 


erg cm -2 s" 1 


erg s _1 




PL 


2.44±0.85 


0.54±0.02 




4.55x10-" 


1.98x10" 


1.21/830 


PL+BB 


7.1 + 1.6 


0.42±0.02 


0.242±0.024 


4.75x10-" 


2.08x10" 


1.02/828 


PL+MEKAL 


9.4±1.5 


0.35±0.08 


1 7 +1 - 9 

1 -'-0.4 


4.79x10-" 


2.11x10" 


1.04/828 


PL+DB 


11.3+2.3 


0.41+0.04 


0.36±0.06 


4.77x10-" 


2.10x10" 


1.02/828 



(1) For definition of spectral models see text. (2) Observed 0. 2-10.0 keV flux. (3) Source intrinsic X-ray luminosity in the 0.2-10.0 keV band 
(corrected for absorption) for a distance to the SMC of 60 kpc (Hilditch et al. 2005|). 
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